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In this paper we study a new family of black hole initial data sets corresponding to distorted "Kerr" black holes with moderate 
rotation parameters, and distorted Schwarzschild black holes with even- and odd-parity radiation. These data sets build on 
the earlier rotating black holes of Bowen and York and the distorted Brill wave plus black hole data sets. We describe the 
construction of this large family of rotating black holes. We present a systematic study of important properties of these data 
sets, such as the size and shape of their apparent horizons, and the maximum amount of radiation that can leave the system 
during evolution. These data sets should be a very useful starting point for studying the evolution of highly dynamical black 
holes and can easily be extended to 3D. 

ner |Q and Brill and Lindquist . At about the same 
time Brill developed an initial data set for gravitational 
wave spacetimes that involved specifying a free function 
in the conformal part of the 3-metric that determined 
the distribution of the gravitational wave energy in the 
spacetime JlGf . These time symmetric initial data sets 
were combined to produce the NCSA "Brill wave plus 
black hole" spacetimes, generalizing the single wormhole 
to a highly distorted black hole spacetime The 
evolution of these data sets has been very useful in un- 
derstanding the dynamics of highly distorted black hole 
spacetimes, and also in preparing for studies of collid- 
ing black holes, as colliding holes form a highly distorted 
single hole immediately after the merging process. Rotat- 
ing, stationary black hole data sets, which are not time- 
symmetric, were first discovered by Kerr p8| in 1963, and 
then a nonstationary form was discovered by Bowen and 
York ll^. The conformally flat Bowen and York con- 
struction for single black holes was then generalized to 
multiple conformally flat holes with angular momentum 
by a number of authors (for example, see [Q) leading to 
recently computed 3D data sets for two black holes with 
arbitrary spin and angular momentum pl[ . 

In this paper we discuss a new family of distorted ro- 
tating black hole initial data sets, some of which were 
evolved and studied in Ref. |2^j2^,Q. Although the 
constructions we present here are axisymmetric, our ap- 
proach is easily carried out in 3D. These data sets com- 
bine a number of the ideas discussed above, generalizing 
the Bowen and York construction by including a "Brill 
wave" so that highly distorted, rotating black holes with 
gravitational radiation can be studied. In addition to 
these new data sets to be described below, this family in- 
cludes as special cases all previous single black hole data 
sets discussed above, including the Schwarzschild, Kerr, 
Bowen and York, and NCSA "Brill wave plus black hole" 
spacetimes. 

As we have shown in Ref. |^,^,^ , these data sets are 
interesting in their own right as dynamical black holes. 
These black holes radiate both polarizations of the grav- 
itational wave held, and their horizon geometries can be 
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I. INTRODUCTION 

There is increasing interest in the study of black 
holes from both astrophysical and theoretical points of 
view. Coalescing black hole binaries are considered an 
important source of gravitational waves for LIGO [Q, 
VIRGO Q and other gravitational wave observatories. 
At the same time, numerical studies of black holes are 
advancing to the point where highly distorted, axisym- 
metric black holes [||J^ and black hole collisions |^j6) can 
be simulated, extracting important physics such as the 
gravitational wave forms emitted and the dynamics and 
properties of the event and apparent horizons [0J|] . Re- 
cent progress has also been made in 3D |^ black hole 
spacetimes. These are some of the reasons it is impor- 
tant to develop a series of initial data sets to evolve. 

However, black hole initial data sets are interesting in 
their own right, apart from evolutions. They are poten- 
tially useful for studying problems related to cosmic cen- 
sorship, such as the Penrose inequality Furthermore, 
it is possible to study greater numbers of data sets with 
higher resolution than is practical to evolve numerically. 
This can lead to important observations of properties of 
apparent horizons and can help to identify particularly 
interesting spacetimes for evolution studies. Initial data 
sets of the type described in this paper can be used to 
analyze the most extreme types of distortion and may 
have bearing on the hoop conjecture [0|2j. The study 
of an apparent horizon in an initial data set may pro- 
vide valuable insight into which spacetimes might have 
the most distorted event horizons. 

Black hole initial data sets, based on the Einstein- 
Rosen bridge construction , have been developed over 
the last few decades. In the early 1960's a series of 
time symmetric, conformally flat wormhole data sets cor- 
responding to two black holes were developed by Mis- 
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so distorted they cannot be embedded in a flat Euclidean 
space. In future papers we will continue to explore their 
dynamics in detail. Moreover, just as in the non-rotating 
case, these data sets can be used to understand the late 
time behavior of coalescing, multiple black hole systems, 
as highly distorted "Kerr" black holes will be formed in 
the process. These data sets can be thought of as repre- 
senting initial conditions for the late stages of that pro- 
cess, and should prove to be a valuable system for study- 
ing it without having to evolve the orbits leading up to 
the merger. 

The paper is organized as follows: In section |l| we 
detail the mathematical construction of the data sets and 
discuss some of their properties, such as their masses, for 
several subclasses of the spacetimes. Then in section [II 
we discuss a series of tools we have developed to study the 
initial data sets, such as the apparent horizons and their 
intrinsic geometry. In section IV we survey many initial 
data sets, describe their features, and discuss application 
of the hoop conjecture and Penrose inequality. Finally in 
section ^ we summarize the results. 



II. CONSTRUCTING DISTORTED, ROTATING 
BLACK HOLES 

A. Initial Data 

In this section we provide details of the construction 
of this new family of distorted, rotating black hole data 
sets. Some details were provided in Ref. |Q, but as this 
paper focuses exclusively on the initial value problem we 
present a more comprehensive treatment here. To moti- 
vate the problem, consider first the Kerr metric, written 
in Boycr-Lindquist form 
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where m is the mass of the Kerr black hole, a is the angu- 
lar momentum parameter, and a superscript denotes 
the Kerr spacetime. 

We would like to put this in a form that is free of coor- 
dinate singularities, and that has properties similar to the 
form used in the non-rotating black hole studies |0,n^,E^ . 



Generalizing the transformation used in Ref. p 17 27|] we 
define a new radial coordinate rj through 

r r+ cosh^ (r]/2) - r_ sinh^ {■q/2) , (2a) 
r± — m :k \/ m? — a^. (2b) 



With this transformation the spatial part of the Kerr 
metric can be written as 

df = [e~2* {df]^ + dO^) -f sin^ ed4>'^] (3) 



where 
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Thus we see that our coordinate transformation (||) has 
resulted in a "quasi-isotropic gauge" |2^ for the Kerr 
spacetime. Notice also that if a = then go = and 
we recover the Schwarzschild 3-metric. This metric is 
now in the form used in the previous NCSA black hole 
studies , although here the functions go and 

are determined by the Kerr spacetime. 

One may check that the 3-metric defined by Eq. (|^) 
is invariant under the transformation rj —rj, i.e., 
there is an isometry operation across the throat of an 
Einstcin-Roscn bridge located at 77 = 0, just as in the 
Schwarzschild spacetime [Q. This construction has two 
geometrically identical sheets connected smoothly at the 
throat, located at 77 = 0. In terms of the more familiar 
radial coordinates, this condition can be expressed as 



/(4r-), 
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where r is a generalization of the Schwarzschild isotropic 
radius, defined as 



(6a) 



r is related to the usual Boyer-Lindquist radial coordinate 
via 



(6b) 
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Note that in the Kerr spacetime the horizon, located at 
r = m + V — a^, is at f = V — a^/2 in the r co- 
ordinates, or at 77 = 0, just as in previous studies of the 
Schwarzschild spacetime pO| ]. The isometry condition 
77 —77 will be used in the construction of more gen- 
eral initial data sets to be described below, and is also 
imposed during the evolution as well Q . 

The 3-metric, however, comprises only half of our ini- 
tial data. The extrinsic curvature on the initial slice is 
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required as well, and for the Kerr spacetime it is given 
by the equations 
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where we note that K., 



^'q Kij, and all extrinsic cur- 



vature components excepting the two listed are zero. 

One can show that this metric and extrinsic curvature 
satisfy the both the Hamiltonian constraint 



R + K^ - K'^K,. 







(8) 



which breaks the stationary nature of the spacetime. As 
we will show below, appropriate choices of this function q 
(along with appropriate solutions to the momentum con- 
straint for the extrinsic curvature terms) can lead also to 
Schwarzschild, the NCSA distorted non-rotating black 
hole (as in Ref . jl^ ) , the Bowcn and York rotating black 
hole [M, or a distorted Bowen and York black hole. 

The function q, representing the Brill wave, can be cho- 
sen somewhat arbitrarily, subject to symmetry conditions 
on the throat, axis and equator, and falloff conditions at 
large radii ■ The function q will be chosen to have 

an inversion symmetric (i.e., symmetric under rj —rj) 
Gaussian part given by: 



and the three momentum constraints 



V, (K'^ ~ -f'^K) = 0, 



(9) 



as they must. 

Finally we must specify the lapse and shift to com- 
plete the system for evolution. The choices that make 
the stationary nature of this Kerr spacetime explicit are: 
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We have now written the Kerr spacetime in 3-1-1 form, 
in a coordinate system like that used to evolve the non- 
rotating distorted black holes in previous work. This is a 
stationary spacetime, and it can be evolved using a dy- 
namic slicing condition just as the Schwarzschild space- 
time has been evolved using maximal and other slicing 
conditions. But just as in the case of Schwarzschild, we 
can distort this rotating black hole with a gravitational 
("Brill") wave and study its response. 

We accomplish this by adding a free function ^(t?, 6) 
to the conformal 3-metric in Eq. (||), just as in Ref. p7| . 
Then the conformal factor ^'o, which is a solution only 
for the stationary Kerr spacetime, will no longer satisfy 
the Hamiltonian constraint (^). Instead we must solve 
this constraint equation for the general conformal factor 
'i' as we discuss below. This generalization results in a 
3-metric with the following form: 



{drf + dO^) 
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Apart from the function go this is the same 3-metric used 
for the non-rotating distorted black hole spacetime. This 
illustrates the fact that unlike the Schwarzschild space- 
time the undistorted Kerr spacetime is not conformally 
flat [R 7^ 0, where R is the scalar curvature of the con- 
formal part of the 3-metric). We may, in fact, make the 
spacetime conformally flat by adding a "Brill wave" ||l^] 
specified by q = go- This distorts the Kerr black hole, 



S± = (??± 770) /cr^. 
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This form of the Brill wave will be characterized by sev- 
eral parameters: Qq (its amplitude), a (its width), 770 
(its coordinate location), and n, specifying its angular 
dependence, which must be positive and even. We note 
that Eqs. ( p^ simply provide a convenient way to pa- 
rameterize the initial data sets, and to allow us to easily 
adjust the "Brill wave" part of the initial data. Many 
other devices are possible. To narrow the field of study, 
in this paper we restrict ourselves to sets characterized 
by n = 2, and rjo = Q for all data sets, and will commonly 
use cr = 1.5 

With this form of the Brill wave function q, a num- 
ber of black hole data sets can be constructed. As dis- 
cussed above, given the appropriate extrinsic curvature 
choice, the stationary Kerr solution results if g = 0. For 
other choices the Bowen and York form results if g = go , 
a distorted Bowen and York spacetime can be made by 
setting q — sm^ Oqo + go, and a distorted Kerr space- 
time can be made by setting q — sin" Oqc- Furthermore, 
for spacetimes without angular momentum this metric 
form contains both the NCSA Brill wave plus black hole 
spacetimes and, as we will see below, a new black hole 
spacetime corresponding to odd-parity distortions of the 
Schwarzschild spacetime. 

As a consequence of this generalization of the metric, 
we must now solve both the Hamiltonian and momentum 
constraints. The Hamiltonian constraint equation (^) can 
be expanded in coordinate form to yield 
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where we have assumed the initial slice is maximal, i.e. 
trK — 0. There are in principle three momentum con- 
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straints to be satisfied by the extrinsic curvature compo- 
nents. However, in our spacetimes we require a "time- 
rotation" symmetry about the initial slice, or invariance 
of the metric under (t, (/)) — + (— — 0). Writing the ex- 
trinsic curvature tensor as 
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this time-rotation symmetry requires that H^ = Hg — 
He = Ho ~ 0, leaving only He and Hp to be deter- 
mined. This also automatically makes the initial slice 
maximal. 

Under these conditions the only nontrivial component 
of the momentum constraints is the ((> component: 



O^He sin" 



9 + dg {Hf sin^ 9^ 



0. 



(16) 



Note that this equation is independent of the functions q 
and go. This enables us to choose the solutions to these 
equations independently of our choice of metric pertur- 
bation. 

At this point we may regard Eqs. ( |l^jl^ ) as defining 
the problem we wish to solve. The boundary conditions 
for these equations, which are extensions of those used in 
Ref. [^ , may be summarized as: 

Isometry Conditions: Across the throat of the black 
hole, labeled by = 0, we can demand a condition that 
the spacetime has the same geometry as ry — ?/. We 
can require this not only for the initial data, but also 
throughout the evolution of the spacetime. This condi- 
tion may be stated in one of two ways (in axisymmetry) 
to allow for different slicing conditions. Each choice must 
result in a symmetric and an anti-symmetric K^/e to 
be consistent with the Kerr solution. Thus, the form of 
the isometry condition must be different for the symmet- 
ric and anti-symmetric lapse spacetimes: rj —>■ — r; for the 
former, {rj, (p) — > {—rj, —0) for the latter. 

If a lapse that is anti-symmetric across the throat is 
desired, the metric elements with a single 77 index are 
anti-symmetric across the throat, while those with zero 
or two indices are symmetric. The extrinsic curvature 
components have the opposite symmetry of their corre- 
sponding metric elements. The shift is anti-symmetric 
across the throat, while all other shifts are symmetric. 

If a symmetric lapse is desired, the metric elements 
with a single 77 index or single (j) index (but not both) 
will be anti-symmetric at the throat and all others will 
be symmetric. Extrinsic curvature components will have 
the same symmetries as their corresponding metric ele- 
ments. The f3'^ and /?* shifts will be anti-symmetric, and 
the shift will be symmetric. With these symmetries 
enforced, the initial data and all subsequent time slices 
will be isometric across the throat. One can verify that 



all Einstein equations respect these symmetries during 

the evolution if they are satisfied initially. 

Symmetry Axis: All metric elements, extrinsic curvature 

components, and shift components with a single 9 index 

are anti-symmetric across the axis. The remainder are 

symmetric. 

Equatorial Plane: At the equator there are two possible 
symmetries, the Kerr symmetry and the "cosmic screw" 
symmetry (see section [I A4| below) . For the Kerr sym- 
metry, 9 TT — 9, all metric components, extrinsic cur- 
vature components or shifts with a single 9 index are 
anti-symmetric. The remainder are symmetric. For the 
cosmic-screw type boundary conditions the symmetry at 
the equator is {9, </)) ^ (tt — 9, —(f)) and those metric 
elements, extrinsic curvature components, or shift com- 
ponents that contain one 9 or one index (but not both) 
are anti-symmetric. The remainder are symmetric. 
Outer Boundary: At the outer boundary a Robin con- 
dition is used for This condition gives the correct 
asymptotic behavior in the conformal factor to order r~^. 
For the metric given in the form (|l2|) "if has the form 

^ e''/2-|_(77i/2)e~''/^-|-... and therefore obeys the differ- 
ential equation + 29^,5' = 26^*/^. The conformal factor 
is always symmetric at the throat, axis, and equator. 

With this description of the problem and boundary 
conditions, we summarize a number of families of rotating 
black hole solutions that are possible in this formulation: 



1. Kerr and Distorted Kerr 



As shown above, Eqs. ( |l2fjl6| ), together with the 
boundary conditions, are satisfied by the stationary Kerr 
spacetime with the transformations (|^^. In Eq. (p^), 
this spacetime corresponds to q — 0, but as we have 
mentioned above, we may add a Brill wave function 
q = sin" 9qG to distort the Kerr spacetime. From the 
standard ADM expressions for the mass and angular mo- 
mentum of the spacetime ]3l|] we have 



Ma 



DM 
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p- = hi " ^'"^ ^^^^^ 

or in terms of the variables defined in this paper, 

Madm ^ f e"^/^ {~d,,^ + 'i'/2)sm9d9, (18a) 
Jo 

J = P4, = - / -^^He sin^ 9d9. (18b) 



If one evaluates these expressions for the undistorted 
Kerr spacetime, (q = 0), one recovers the expected re- 
sults Madm = fTi (which should be evaluated at 77 = 00) 
and J = am (which may be evaluated at any radius in 
axisymmetry). In practice, we evaluate m at the outer 
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edge of our grid, and J is simply an input parameter. We 
have evaluated and verified J at various radii, as reported 
in previous work |22|. 



2. Distorted Bowen and York 

Perhaps the simplest nontrivial choice of solution for 
the momentum constraint is given by the Bowen and 
York lUl solution 



He = 3 J 
Hp = 0, 



(19a) 
(19b) 



where J is the total angular momentum of the spacetime. 
This solution clearly satisfies the momentum constraint. 
In the language of this paper, the original Bowen and 
York solution has a Brill wave function determined by 



q = qo, 



(20) 



i.e., it is conformally flat, but we may distort it by setting 



qo + sin" Oqc- 



(21) 



This is the solution we often use in constructing distorted 
rotating black hole spacetimes, and a number of them 
have been evolved and discussed in Ref. [|3|. It is some- 
what simpler to construct than the distorted Kerr be- 
cause we can choose our total angular momentum and 
distortion parameters without needing to know the ra- 
tio J/rri^ (an additional value that we would need to 
initialize the Kerr momentum solution). Knowing this 
parameter would require us to know the value to, and 
we will not know that parameter until after we solve the 
Hamiltonian constraint. It is possible, in principle, to 
iterate between guesses for the parameter J jw? and so- 
lutions for the Hamiltonian constraint until arriving at a 
correctly initialized Kerr momentum solution. 

We noted above that the Kerr solution is not confor- 
mally flat, while the Bowen and York solution is con- 
formally flat. Therefore one can consider adding a Brill 
wave to the Kerr solution, making it a conformally flat 
Bowen and York solution. One can study the form of 
the wave functions q and go to understand what sort of 
Brill wave is required to "flatten" the Kerr solution. By 
equating the series expansions of q^ and our Brill wave 
function (c.f. Eq. (|4b|)) 



gapprox = ai sm de 



2 / 2 



a2 sm 



2 / 2 



(22) 



to second order in rj and 9, one obtains relationships be- 
tween the wave shape parameters ai , 02 , cri , and (T2 and 
the Kerr parameters a and m. It turns out that for 
a/m < .9 we find that the sm^ 9 term dominates the 
sin^ 9 term, and cri varies in value from 1.26 to 1.7 for 
rotation parameters 0.01 < a/m < 0.8. Furthermore the 
amplitude ai varies smoothly from zero to about 0.38 for 



all rotation parameters < a/m < 1. We conclude that 
the Bowen and York distortion is similar to a low am- 
plitude n = 2, (T = 1.5 Brill wave distortion on a Kerr 
black hole, which is quite similar in shape and size to the 
distorted, non-rotating black holes evolved previously 

To get a feel for the difference between the qo function 
and its approximation, qapprox, we calculated the ADM 
value of q/Madm as we would for the Kerr spacetime, 
except that we used qapprox instead of qo in the metric. 
We found that the error in a/MADM was less than 1%, 
even for a/m = 0.9. 



3. Odd-parity Distorted Schwarzschild 

Odd-parity gravitational modes, as defined by Regge 
and Wheeler | ]3^ , are metric perturbations with parity 
(—1)^+^ where £ is the angular index of the tensor spher- 
ical harmonic describing the wave. These modes, in ax- 
isymmetry, are formed only from perturbations in 
and gg^. The evolution of these modes does not cou- 
ple linearly to the even-parity modes, and if there is no 
odd-parity disturbance on the initial slice of the space- 
time these modes are never excited even when the full 
non-linear equations are used. 

We can create odd-parity distorted non-rotating space- 
times simply by finding solutions to the momentum con- 
straint equation (|l6|) which drive the evolution of ge^, 
which satisfy the boundary conditions, and yield J = 
when the integral given in Eq. (18b) is performed. 

As a framework for momentum solutions, we consider 
the simple class defined by: 



He = im + f2{v) (4 cos 9 h{9) + sin 9e /3(e)) (23a) 
Hf = -dj2{ri)sin'df3iO). (23b) 

This solves the momentum constraint where /i , /2 and /a 
are arbitrary functions which satisfy the following sym- 
metries: 



hi9) = h{-9) = Mt: - 9) 

f2{v) = /2(-r?) 
lim /2(ry) = 



(24a) 
(24b) 
(24c) 
(24d) 



The Kerr solution is not a member of this class, but the 
Bowen and York solution does fit within it (/i = 3J, 

/2 - /3 = 0). 

Note that any spacetime with fi — will have zero 
angular momentum, as He falls off for large values of the 
radial coordinate 7 7, an d furthermore the angular part of 
the integral in Eq.( |l8b| )is identically zero by construction. 
We require /2 to go to zero only so that we can make the 
spacetime asymptotically Kerr. 

We now turn to a solution that has no net angu- 
lar momentum, but one that is not time symmetric 
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about the initial slice. This solution corresponds to a 
"Schwarzschild" black hole with both even- and odd- 
parity distortions. This generalizes the NCSA Brill wave 
plus black hole spacetimes ||2^ , which allowed only even- 
parity distortions. 

For this we choose the extrinsic curvature components 
such that fi = 0, /2 = (/q, and = 0036* sin" ^^9. 
In order to obey symmetry and regularity conditions, n' 
must be odd and have a value of at least 3 (This was 
arranged to be suggestive of a radiation ^-mode. Using 
n = 3 gives almost pure ^ = 3, while using n = 5 gives 
almost a pure i = 5 component [^.). 

As this data set contains the odd-parity polarization of 
the gravitational wave field, but no angular momentum, 
it is simpler to evolve than than the rotating black holes 
described above, but has richer physical content than the 
standard Brill wave plus black hole data sets. A few 
of these data sets have been evolved and discussed in 
Ref. ||2^, where both even- and odd-parity waveforms 
were extracted. 



4- Cosmic Screws 

There is another solution related to the odd-parity dis- 
torted Schwarzschild data set discussed above. While 
that data set has equatorial plane symmetry, we can al- 
ternatively construct a solution that has something like 
an opposite sense of rotation below the equatorial plane. 
As discussed above, this symmetry is {9, 0) (tt— 0, —(f>), 
thus requiring He to be antisymmetric and Hp to be 
symmetric across the equator. 

Axisymmetric data sets describing the collision of two 
counter-rotating black holes have been dubbed "cosmic 
screws." As our data sets contain only one black hole, it 
can be considered to represent a later stage of evolution 
of cosmic screws, where the holes have already merged. 
Such data sets should be useful in studying the late time 
behavior of colliding spinning black holes, and will be 
evolved in a future paper. The formalism of the previous 
section carries over, but now the symmetry conditions 
on the equator are reversed: fi{9) — — /i(7r — 9) and 

f3iO)=M7T-9). 

With these considerations in mind, we see that /i = 
3Jcos6', and /2 = /s = is suggestive of a cosmic screw 
spacetime (in this case we do not claim that J is an an- 
gular momentum, but use this notation because it is sug- 
gestive of the angular momentum of one of the colliding 
black holes in a cosmic-screw spacetime). 

We can also choose /i = 0, f2 = Qq, and /a = 
sin" 9 which is similar to our odd-parity distortion of 
Schwarzschild except that the equatorial boundary con- 
ditions are like those of the cosmic screw. Here n' should 
be even and have a value of at least 2 in order to obey 
symmetry and regularity requirements. 



B. Numerical Techniques 

To solve the Hamiltonian constraint equation we use 
a Newton-Raphson iteration procedure. The conformal 
factor in the constraint equation is linearized about a 
trial function. The initial trial function, ^q, is usually 
chosen to be 2 cosh(77/2), because that is the solution for 
a spherical black hole with mass AI — 2. At a given 
iteration i we assume the solution has the form ^ = 
^'i + S'i', where S'^ is a small correction function. We 
then evaluate the residual of this guess for '^i, 



{d?,{q-qo) + d'g(q~qo)) 



Hi sin^ 9 + Hj. 



(25a) 



and use it as a source term for the correction function 



AS^ + — (9,2 {q - qo) + dl (q - qo)) 



*; 



Hi sin^ I 



Hi - -e, 



(25b) 



Eq. (25a) is used merely to evaluate e^, but Eq. ( p5b| ) 
is solved with a numerical elliptic solver for S'i'. The 
boundary conditions on the grid are given as follows: 6"^ 
is symmetric across the throat, the axis, and the equa- 
tor. At the outer boundary we use the Robin condition, 
djfS'i' — —^S'i>. Typically we use e ~ 10^^ as our stop- 
ping criterion. 

Once (5* is obtained, we try again, with ^'i+i = + 
(5*. We repeat this procedure until the a solution to Eq. 
( p^ ) with appropriately small error (maximum value of 
ei on the grid) is achieved |Q . 

Some of the data sets considered in this paper use Brill 
waves that are very wide. Because of this, one may worry 
whether the outer boundary of the grid is sufficiently 
large to get an appropriate asymptotic behavior for the 
Hamiltonian constraint. To determine whether the posi- 
tion of r]max is valid, we measure the quantity 



^ = 9,, log 



5- 



(26) 



which should approach — 1 /2 in a spacetime that asymp- 
totically approaches the Kerr spacetime. We adopt the 
criterion that ^ must be equal to this value to within 
0.1%. If it is not, we regard our choice for value of r]max, 
the outer boundary of our grid, to be too small and adjust 
accordingly. 



C. Limits on tlie magnitude of Brill wave amplitudes 

Cantor and Brill ||3^ showed that for a certain class of 
spacetimes there are limits on the value of the amplitude 
of the Brill wave in a spacetime with topology R^. Their 
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basic technique can be applied to these rotating black 
hole spacetimes, which have topology RxS^. We present 
a simplified version of this analysis. 

The strategy is to derive an inequality that contains 
an arbitrary function /, but docs not explicitly contain 
the conformal factor '5, and which requires asymptotic 
flatness of the spacetime as one of its conditions for va- 
lidity. We use the freedom in / to locate Brill wave pa- 
rameter spaces which fail to satisfy the inequality. When 
the inequality fails, it is then clear that the assumption 
of asymptotic flatness has been violated, and the value 
of the input Brill wave parameters which specify it are 
outside our domain of interest. 

For this section, we define the metric as 



U4/ 



Brill started by considering the volume integral 



V/- V/dT/ 



(27) 



(28) 



where / = u^*, and where m or / is an arbitrarily chosen 
function. We can expand the above integral for / in 
terms of u and VP: 

V • {u^^V^^ - u^^A^ = 2u^Vu ■ + u^V^ ■ V^. 

(29) 

We have now obtained both a surface and a volume piece 
for the integral. Because we want a final equation which 
does not include 4*, we discard the terms under the in- 
tegral which still contain it explicitly after our transfor- 
mation. Since these terms are manifestly positive, the 
remainder must be less than our original integral. We 
now have an inequality. Next we require that ^ not van- 
ish (this is only necessary if we have non-zero extrinsic 
curvature), is finite, and that Vu go to zero faster than 

To evaluate the surface terms we make use of the 
asymptotic properties of VP, namely that ^ = 1 + ^ + 
0{r~^) at large radii. The surface term is 



M^^-V* • dS. 



(30) 



This last expression will be —2n times the ADM mass 
of the system if we choose u '^^^ sX, rj ^ oo. If 
we choose m — > as 77 — > cx) then the integral will not 
evaluate to the ADM mass (in fact it will obtain zero). 
When we perform this surface integral at the throat of 
our spacetime we exploit the isometry 



= 0, 



(31) 



to obtain the final form for our inequality 



"^T^fn^ooMADM 



\7f-<7f+-fR]dV> 



TT / f^^Qsin9 d9. 
'0 



(32) 



where we have assumed that / is independent of 9 at 
large 77. 

This inequality can only fail if our assumptions about 
the form of ^E" are incorrect. Therefore, when it does 
fail, we can conclude that an asymptotically flat metric 
cannot be formed by the given Brill wave parameters. 

We explored the regime where q = sin" 9 qc + Qo = 
2Qq sin^ 9 e~'^ /'^ and found two forms for / which pro- 
vided us with good bounds on the integral. These forms 
are 



and 



/ = exp(-e77^ - ^77/2) 



/ = cos'' 9 exp {—erj^ — M^/2) 



(33) 



(34) 



Replacing the inequality with an equality in Eq. (p^), 
performing the integral, and then solving for Qq provides 
us with a function, QQ{^,e,a), whose range supplies the 
forbidden amplitude values Qq and whose domain is e > 
0, cr > 0, and is a real number. 

The functions Qo{n, e, a) that we obtained were able to 
provide both an upper and a lower bound for the allowed 
values of the input parameter Qo because they trace out 
a function similar to an hyperbola in fi space for given 
values of a and e. The behavior of one function is il- 
lustrated in Fig. 1^. In this figure we see two distinct 
extrema. One extremum is a greatest possible negative 
value for Qo and provides us with a lower bound for this 
parameter. The other is a minimum in the possible pos- 
itive values for Qq and this provides us with a maximum 
value for Qo- 



Finding Min/Max 




FIG. 1. We plot the range of the Brill wave amplitude func- 
tion Qowhich are inconsistent with asymptotically flat solu- 
tions of the Hamiltonian constraint for a given "trial function" 
(see text). In this example plot, the values outside the range 
-2.317 < Qo < 56.24 are excluded. 
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In Table | we have give the strictest upper bound for 
the function Qq that we were able to find using the above 
inequality, Eq. (^), and trial functions, Eqs. (^,|4|). 
These bounds, together with the parameters /i and e 
which were chosen to minimize them, are given in the 
table. 



(J 


e 




n 


Oo 


0.25 


1.726 


-2.928 





3.283 


0.50 


1.078 


-3.186 





5.122 


0.75 


0.8455 


-3.64 





7.516 


1.00 


0.7042 


-4.037 





10.74 


1.25 


0.6171 


1.096 


4 


11.97 


1.50 


0.5580 


1.097 


4 


10.99 


1.75 


0.5049 


1.087 


4 


10.35 


2.00 


0.4587 


1.076 


4 


9.914 


3.00 


0.3300 


1.046 


4 


8.993 


100.0 


0.0100 


1.014 


4 


7.6113 



TABLE I. Upper limits on Qo- The parameters e, /i, and 
n are arbitrary parameters were chosen in an attempt to find 
the smallest upper bound for the allowed values of Qo, the 
wave amplitude, as a function of the wave width, a. 

If one explores the n ~ I form of this distortion and 
takes the limits tr — *■ oo and e — + in that order one 
arrives at the amplitude Qo = 7.577. Thus, it seems that 
while there are more strict limits on larger width waves 
in general, there is no reason to suspect that any width 
is too large. 

In Table || we have the strictest lower bound for the 
function Qo that we were able to find using the methods 
of this section. 



a 


e 




n 


Qo 


0.25 


0.0025 


2.81 





-1.261 


0.50 


0.0025 


2.31 





-1.437 


0.75 


0.0558 


1.94 





-1.486 


1.00 


0.1053 


1.62 





-1.456 


1.25 


0.1133 


1.44 





-1.404 


1.50 


0.1081 


1.33 





-1.351 


1.75 


0.1001 


1.26 





-1.303 


2.00 


0.0921 


1.21 





-1.260 


100.0 


0.0487 


1.00 





-0.765 



TABLE II. Lower limits on Qo. The parameters e, n, and 
n are arbitrary parameters were chosen in an attempt to find 
the greatest lower bound for the allowed values of the wave 
amplitude Qo, and the wave width, a. 

It is interesting to note, however, that the pure Bowen 
and York spacetime (which is parameterized by J, not Qo 
and can be regarded as a distortion of the Kerr space- 
time) can be specified apparently without limit. This 
is consistent with the calculation done here, because the 
conformal Ricci scalar vanishes for conformally flat space- 
times such as the Bowen and York spacetime. This as- 



sures us that the inequality given in Eq. (|3^) is trivially 
satisfied. Numerical evaluations of the Bowen and York 
black hole (see Q]) seem to get asymptotically closer 
to the extremal Kerr spacetime as J is increased, i.e. 
J/m^ = a/m gets ever closer to unity confirming this 
observation. 

We note also that we can use Eq. ( ^^ to obtain a lower 
bound on the ADM mass for numerically derived data 
sets. Using / = 1 we obtain 



1 r 

Madm > 1 + t / (g - go) sin 6* dO 
4 Jo 



(35) 



where the integral is evaluated on the throat. We calcu- 
late this after computing 5* on the initial slice to deter- 
mine whether the ADM mass of the spacetime is consis- 
tent with asymptotic flatness. 



III. APPARENT HORIZONS 

A. The numerical method 

We will use the same method for finding the AH (ap- 
parent horizon) as described in previous work Jl0| , p3t , 
where surfaces that satisfy the trapped surface equation 



D,s' + K,,s's^ -K = 



(36) 



are sought. In the above, s* is an outward pointing unit 
normal vector which defines the surface. For initial data 
sets of the type considered in this paper, the terms con- 
taining the extrinsic curvature drop out, and therefore 
finding the apparent horizon reduces to finding an ex- 
tremal surface. We will, however, be using a different 
strategy for searching than we used when we were evolv- 
ing our spacetimes p^ , p3t . For extremely distorted data 
sets we found that multiple trapped surfaces may ex- 
ist, and it is important to find the outermost surface. 
Therefore, we use a "brute force" approach and repeat- 
edly call the solver, using many constant radial zones 
covering the inner portion of the grid as an initial guess. 
After the horizon finder has finished executing for each 
initial guess, the solution with the largest radial coordi- 
nate position that satisfied the tolerance of the solver was 
determined to be the apparent horizon. 

This strategy relies on our previous experience, which 
has shown us that even horizons with highly distorted 
geometry are nearly spherical in coordinate space, in our 
choice of coordinates. Moreover, we know, from Gib- 



bons 1 35 1 , that our apparent horizon must have the topol- 
ogy of a 2-sphere. As in previous papers H^,^, we will 
study the intrinsic geometry of the black hole horizon as 
a measure of its distortion. 



B. Measures of the horizon 

There are a number of geometric quantities that can 
be used to measure the apparent horizon that provide 
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insight to the physical characteristics of the black hole. 
Examples are its shape, area, and local curvature. One 
straightforward measure of a black hole's shape in ax- 
isymmetry is the ratio of its polar circumference, Cp, to 
its equatorial circumference, Cg. This ratio has been used 
previously, but we point out that it has an inadequacy. 
If used to study the surface of an event horizon, this 
number may be infinitely large for some configurations; 
one need merely consider the case where two equal mass 
black holes are colliding and are just touching so that 
Ce vanishes. Although we do not expect this problem 
to affect the apparent horizon, however, we would like 
to use the same measure of oblateness for both apparent 
and event horizons, one which better captures the intu- 
itive concept of oblateness. One such measure is Cp/C^, 
where Ca = 2tt{2^J A/lQir) is the "areal circumference." 
Note that the quantity Cp/CA tends to be much smaller 
than Cp/Ce- For the case of a c = 3.5, Qo = 5.64, 
black hole apparent horiozon we have Cp/CA — 85 and 
Cp/Ce = 19 X 10^. 

Another important quantity that can be obtained from 
the horizon is its mass 



a set of black holes with a 
amplitude, Qq, in the range - 



= 1.5, and the Brill wave 
-1< Qo < 7. 



A 

16^ 



(37) 



and the related quantity, the MRL (maximum radiation 
loss), 



MRL = 1 - Mah/^'I 



(38) 



The MRL is the amount of energy that can be emitted 
if the black hole swallows none of the radiation energy 
present outside the horizon on the initial slice during its 
subsequent evolution. 

We will also be interested in discussing what we call 
the "transition points" of the horizon. These are points 
in parameter space of the initial data problem at which 
the apparent horizon ceases to be on the throat (which is 
always a trapped surface) of the spacetime. The transi- 
tion points can be accompanied by discontinuous changes 
in the character of the horizon, as measured by the quan- 
tities discussed here. 



IV. SURVEY OF INITIAL DATA SETS 

We now apply a number of these horizon measures de- 
veloped above to analyze a series of initial data sets, be- 
ginning with distorted, non-rotating holes. 

A. Time-symmetric Brill wave data 

First we begin by examining the shape parameter, 
Cp/Cyi, for a series of Brill wave plus black hole space- 
times. These non-rotating black holes provide a good 
starting point for our survey. Fig. |^ shows Cp/CA for 
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FIG. 2. We plot the apparent horizon shape parameter 
Cp/Ca, where Cp and Ca are the polar and "areal" circumfer- 
ences of the surface, for a set of non-rotating distorted black 
hole initial data sets. The shape parameter is plotted against 
Qo, the amplitude of the distorting Brill wave. All spacetimes 
shown here have the same Brill wave width, a = 1.5, the width 
that best characterizes the geometry of the Kerr spacetime. 
This plot shows the same data as computed on three grids 
with increasingly more zones, which illustrates how close we 
are to the exact solution. These results are second order con- 
vergent. 

The apparent horizon generally becomes prolate if it is 
distorted with a positive amplitude Brill wave, and oblate 
if it is distorted by a negative amplitude Brill wave . 
However, the shape of the horizon undergoes two sharp 
transitions, one at about Qo = —.65 ||l^ and the other 
at about Qo — 4.5. Between these two values of Qq 
the horizon is on the throat, and the ratio Cp/CA grows 
exponentially. The further we move outside this region 
the further out the apparent horizon moves, and the more 
spherical the black hole becomes. 

The spacetime behaves as if increasing amounts of en- 
ergy were being deposited in the region somewhat away 
from the horizon as jQol is increased. Beyond a certain 
critical threshold, one would expect a new horizon would 
be formed. Beyond this range, the spacetime is still be- 
coming more distorted, but the distortion is driving the 
apparent horizon further out on the grid and as a result 
the distortion is now primarily interior to the apparent 
horizon. Thus the exterior region becomes more spherical 
as this parameter is increased. 

Next we show the effect of changing the width of the 
Brill wave. In Fig. ^ we show a set of spacetimes like 
those shown above, but a larger width Brill wave is used 
(cr = 3.5). We see that the horizon stays on the throat 
for larger values of Qo, thereby making the maximum 
distortion parameter much larger. This results in a value 
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of Cp/Ca that is about an order of magnitude larger 
than the cr = 1.5 distortion. The lower transition is still 
at about Qo = —-65, but now the upper transition has 
moved to Qo = 5-65. The decrease in Cp/CA as we move 
away to values of Qo larger than 5.65 also becomes less 
sudden. 



Brill Wave Spacetlmes 
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FIG. 3. This plot illustrates the basic properties of the ap- 
parent horizon shape parameter, Cp/Ca, for the large width 
Brill wave, a = 3.5. The apparent horizon becomes extremely 
distorted near Qo = 5.65, with C^/Ca ~ 90. This plot shows 
the same data as computed on three grids with increasingly 
more zones. 

In Fig. m we now show the data from the previous two 
graphs along with a number of spacetimes with varying 
Brill wave widths, cr. The values of a considered are 0.5, 
1.0, 1.5, 2.0, 2.5, 3.0 and 3.5 We wiU call this group of 
spacetimes set la. We see that the function Cp/CA of Qo 
changes only slightly as a is varied for a large portion of 
the plot. There is a sharp transition that occurs for pos- 
itive amplitude distortions which is strongly dependent 
on the width of the wave. Beyond this transition point 
the AH ceases to be on the throat. There is another tran- 
sition point for the negative amplitude Brill wave space- 
times, and for all but the narrowest of the Brill waves 
considered (cr = 0.5) it lies at about Qo — —0.65. The 
region between these two transitions constitutes what we 
will call the "main branch" and is relatively independent 
of cr. The "main branch" for the range 1.8 < Qo < 5.6 is 
approximated to an accuracy of 1% or less by 

Cp/Ca = exp (-0.287055 + 0.728645 Qo + 0.019728 Qo^) , 

(39) 

and is plotted by a heavy dashed line in Fig. |4[ Since 
the horizon is simply the throat (77 — 0) for the "main 
branch" spacetimes it is not surprising that log{Cp/CA) 
should depend linearly on Qo as it is simply reflecting 
the effect of Qo upon the conformal metric. 



Main Branch 




FIG. 4. We illustrate the prolateness (large Cp/Ca) of the 
apparent horizon for spacetimes with Brill wave distortions 
of varying amplitude, Qo, and a variety of widths, cr. Also 
plotted is a fit of the "main branch," a curve on which horizons 
lie for all distorted black hole spacetimes considered. For 
each family we consider, there is a transition point, a critical 
amplitude above which the horizon jumps off the throat and 
becomes more spherical. 

Also worthy of note is that above cr = 2.0 the upper 
transition point becomes nearly fixed at Qg ~ 5.65. Some 
small variation occurs, but not much. Roughly cr = 2 is a 
width above which the most prolate spacetimc occurs at 
about Qo = 5.6 with a shape characterized by Cp/CA ~ 
88. Likewise, it appears that cr = 1 is a width above 
which the most oblate spacetime occurs at about Qo = 
— .65 with a shape characterized by Cp/CA ~ -77. 

In Fig. H we see the effect of Qo and cr on the MRL 
(maximum radiation loss) of the spacetimc for set la. 
For reference we have marked the positive Qo transition 
points, the peaks in the plot of Cp/CA that are seen in 
Fig. |, with diamonds in Fig. |. The MRL mcreases, m 
all cases, as we increase |Qo| until it reaches some max- 
imum value. As we increase |Qo| further, we arrive at a 
transition point. Along the path of increasing |Qo|, the 
energy of the spacetime first increases primarily because 
we are adding wave energy outside the horizon. As we 
travel further in this parameter space the mass increases 
primarily because energy is being added to the hole itself. 
Another feature that can be observed in Fig. |5| is that the 
MRL of the transition point begins to rise much more 
rapidly with increasing Qo after the Brill wave width be- 
comes greater than 2. 

In both Fig. ^ and Fig. || we have plotted the curves for 
Brill waves with widths cr = 1.5 and a = 3.5 respectively 
at three different grid resolutions: 600 x 96, 300 x 48, 
and 150 x 24. We note that in both cases the resolution 
lines are very close together and that the convergence 
exponent, ^, as defined by, 

{7^) -{7^) + - (40) 

\^A/num \^ A J true 
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is very nearly 2 at all points on both graphs. The larger 
width Brill waves were more difficult to calculate, but our 
code is capable of calculating these spacetimes accurately. 



Effect of Brill wave width on M RL 
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FIG. 5. This plot illustrates the maximum radiation loss, 
or MRL, for spacetimes with Brill wave distortions of varying 
amplitude, Qo, and a variety of widths, a. The diamonds 
mark the location of the maximally prolate horizon for each 
curve 
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FIG. 6. This plot illustrates the behavior of the negative 
amplitude, narrower width Brill wave spacetimes. The be- 
havior here is qualitatively similar to the larger width posi- 
tive amplitude spacetimes depicted in Fig. 4, except that the 
apparent horizons are all oblate. As Qo becomes more nega- 
tive the apparent horizon jumps off the throat at a transition 
point and the spacetime becomes more spherical. 



B. Negative amplitude Brill wave distortions 

In this section we focus on the region of the time- 
symmetric distorted black hole spacetime in which the 
amplitude of the distortion is negative. In this regime we 
find a behavior that is qualitatively similar to the behav- 
ior seen for positive amplitude distortions in set la, but 
we found that we needed to use narrower Brill waves to 
explore this effect for negative amplitude spacetimes. If 
we consider Fig. ^ we see that all but the cr = 0.5 curves 
are very similar for the negative amplitude portion of the 
graph. Because of this we use widths cr = 0.3, 0.4, 0.5, 
... 0.9 and we plot this data in Fig. ^. We label this data 
as set lb. 

Note in Fig. ^ that the maximum distortion (mini- 
mum Cp/Ca) is increased by widening the Brill wave 
until a limiting value is reached, a behavior similar to 
that observed for positive amplitudes in Fig. |4[ There 
is also a continuation of the main branch for the seg- 
ment of the spacetimes in which the wave is on the 
throat, although we have not computed a fit in this case. 



C. Brill wave distortions of Bowen and York 
rotating black holes 

In this section we discuss spacetimes with a — 1.5, 
— 1 < Qo < 7 and angular momentum parameter J — 
0.25, J = 0.5, and J — 5, and we will call this collection 
set 2a. We plot Cp/CA for this set in Fig. |^. One can 
see from the figure that the rotating black holes still lie 
along the main branch, but the transition points occur 
earlier (for smaller Qo)- 

Adding angular momentum to Brill wave spacetimes 
reduces the value of Cp/CA, thereby making the space- 
time more oblate. For the larger amplitude Brill waves 
this effect is quite marked. In Fig. || we make this point 
clearer by graphically depicting the embeddings of sev- 
eral distorted black hole spacetimes. Three Brill wave 
spacetimes are taken as starting points: Qo = —0.1, 0.0 
and 0.1, and are plotted with thick lines. Next to each 
of these lines is a thinner one depicting the shape of the 
horizon after setting J = 10. This general trend, that 
angular momentum makes an horizon more oblate was 
a feature that was observed previously by Smarr for the 
Kerr spacetime p^ . 
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Effect of Rotation 
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FIG. 7. This figure illustrates the effect on the apparent 
horizon of a Brill wave distortion on a Bowen and York black 
hole. The addition of angular momentum makes the tran- 
sition point, the point where the horizon leaves the throat, 
occur sooner, although for small values of Qo the value of 
Cp/Ca is very close to the "main branch." 
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FIG. 8. In this plot we see the embeddings of several dis- 
torted black hole spacetimes. The heavy lines are spacetimes 
distorted by Brill waves: a prolate spacetime with Qo ~ 0.1, 
an oblate spacetime with Qo = —0.1, and a Schwarzschild 
spacetime with Qo ~ 0.0. The lighter lines represent these 
same spacetimes with angular momentum added. 

As mentioned above, one generally expects that rotat- 
ing black holes will be potentially more oblate than non- 
rotating spacetimes. However, there is an exception to 
this. In Fig. ^we extend the above plots to show in more 
detail the effect of rotation on the negative amplitude 
spacetimes. We see data sets which cover input param- 
eters in the range —0.657 < Qo < —0.655, with cr = 1.5 
and < J < 10 {set 2b). Thus, while the addition of 
angular momentum to a spacetime usually makes an AH 
more oblate, this is not always true. In fact, the addition 



of the Bowen and York distortion makes the minimum 
Cp/Ca as a function of Qo larger. In other words, while 
adding a negative amplitude Brill wave to a given black 
hole can make it more oblate, increasing J also makes 
the horizon jump off the throat at a smaller value of \Qq\ 
so that it does not achieve the maximum oblateness seen 
with no rotation. 



Effectof J on Minimum C /C„ 
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FIG. 9. In this plot we see data sets which cover input 
parameters in the range —0.657 < Qo < —0.655, with a = 1.5 
and < J < 10 [set 2b). It illustrates that although J usually 
makes an horizon more oblate, it can also cause the horizon 
to leave the throat for sufficiently negative values of Qo and 
thereby make it more spherical rather than more oblate. 

Although we do not provide a plot, we note that as the 
parameter J is increased for these spacetimes, the MRL 
decreases. 



D. Brill wave distortions of end-stage "Cosmic 
Screws" 

We next examine a collection of "cosmic screw" space- 
times, referred to as set 3. The data in set 3 is identical 
to the data in set 2a (distorted Bowen and York black 
holes) with the exception that it uses the cosmic screw 
form of the perturbation. The value of Jscrew clearly 
creates a much smaller distortion of the horizon than the 
Bowen and York angular momentum parameter of the 
same value, as is evident in Fig. |l^. Indeed, it is almost 
impossible to distinguish the Jscrew = 0.5,0.25 space- 
times. This is reasonable since this distortion can be 
thought of as two counter-rotating and therefore counter- 
acting Bowen and York distortions. 
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Cosmic Screw Black Hole Horizons 
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FIG. 10. This figure illustrates the effect on the apparent 
horizon of a Brill wave distortion on a black hole with a "Cos- 
mic Screw" distortion (counter-rotating above and below the 
equatorial plane but without net angular momentum, see text 
for details) . As in the case of the Bowen and York black holes, 
the boundary conditions the Cosmic Screw has the effect of 
smoothing out the Cp/Ca curve, although it has a lesser ef- 
fect than a Bowen and York distortion with the a distortion 
parameter of equal value. 



E. Odd-parity distortions of Schwarzschild 

In this section we examine set 4, which comprises odd- 
parity distortions of the form /i = 0, /2 = and /a = 
COS0 (see Eq. (23)), with parameters 770 = 0, < Qo < 
100, cr = 0.5, 1.5, 2.0, 2.5, 3.0, and 3.5. 



This family of spacetimes, depicted in Fig. [Til, displays 
the same property of the pure Bowen and York data dis- 
tortions (in which J, the angular momentum parameter, 
can be thought of as representing the amplitude of the 
extrinsic curvature distortion), namely that there is no 
limit to the magnitude of the parameter, J, that is possi- 
ble, either theoreticaUy or in any numerical test done to 
date. We do not feature the negative amplitude distor- 
tions of these sets because the relevant terms are squared 
in the Hamiltonian constraint, making their overaU sign 
unimportant. 



FIG. 11. This plot shows the prolateness of the horizon, 
Cp/Ca, for a odd-parity data sets with a variety of distortion 
widths, (J. As was the case for the Brill wave distorted space- 
times (compare to Fig. 4 above), the maximum in the Cp/Ca 
curve marks the place where the apparent horizon moves away 
from the throat. 

Unlike the Bowen and York spacetime, for which the 
horizon remains on the throat for very large J, the hori- 
zon moves away from the throat in the odd-parity dis- 
tortion of Schwarzschild for sufficiently large Qo- In this 
respect these data sets are similar to set la. They also 
resemble set ia in that larger width distortions make it 
possible to achieve more distorted horizons. Set 4 is un- 
like set la in that the distortion never becomes even re- 
motely as large (note that this distortion is in the ex- 
trinsic curvature and affects the 3-metric only through 
the conformal factor, it is not the standard Brill wave 
distortion) . 

Fig. |l^ shows the behavior of the MRL function for set 
5. As in set la this data set has its peak MRL somewhat 
before the horizon departs from the throat. It is other- 
wise qualitatively similar to the MRL curves depicted in 
Fig. i. 
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Odd-Parity Schwarzschild MRL 
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FIG. 12. This plot shows the MRL for a variety of 
odd-parity data sets. Quahtatively it has the same features 
observed in Fig. 5, except that the magnitude of the MRL is 
much smaller. As in Fig. 5, diamonds mark the amplitudes 
at which a given black hole spacetime has its maximum value 
of Cp/Ca. 



F. Hoop Conjecture 

An important question that wc wish to discuss is 
"How distorted can a black hole be?" The hoop con- 
states that an event horizon will form if 



jecture |11 



and only if the matter is sufficiently compact in all di- 
rections. This conjecture would seem to suggest that the 
black hole horizon should also be sufficiently compact in 
all directions. Although interesting work has been done 
on the hoop conjecture for black hole spacetimes (see, 
e.g. p7| , ^ ) much remains to be done. 

As we have seen above, black hole apparent horizons 
can be extremely distorted, nonspherical objects as mea- 
sured by the shape parameter Cp/CA- We have shown 
cases where the AH is extremely prolate {Cp/CA ~ 88) 
or oblate {Cp/CA ~ 0.77). However, these results are 
for the apparent horizon, and not the event horizon. A 
careful study of distorted event horizons, and implica- 
tions on the hoop conjecture will be presented elsewhere. 
In this section we restrict our attention to the apparent 
horizons we have found, and discuss whether there is a 
limit to how distorted they can be. 

The fact that there are limits on the amplitude parame- 
ter (Qo) is somewhat suggestive of a limit on the degree of 
distortion possible for the apparent/event horizon of the 
black holes considered here (since the amplitude param- 
eter controls the degree of distortion of the black hole). 
Indeed, it allows one to demonstrate the limit numeri- 
cally for a Brill wave of a given shape (as given by cr, 770, 
and n). 



However, there may also be a good reason to suppose 
that there is no limit on how large or small the quantity 
Cp/Ca may become for an asymptotically flat spacetime. 
We might be able to engineer spacetimes with two grav- 
ity wave distortions: one widely extended, the other more 
comparable to the size of the hole. In effect, the larger 
distortion can be used to "fool" the black hole into think- 
ing that the spacetime is not asymptotically flat. Thus, 
it may be that we can find spacetimes with arbitrarily 
large Cp/CA if we combine distortions on a variety of 
scales. This possibility is not investigated here. 



G. McLximum Prolateness 

It is difficult to pin down a maximum prolateness for all 
the types of spacetimes considered here. It seems reason- 
able, however, to simply focus on Brill wave spacetimes 
since the odd-parity and momentum distortions usually 
make the spacetime more oblate or have small effect. For 
any given value of cr it is a straightforward matter to ob- 
tain the value of Qq which maximizes Cp/CA because 
we can search the entire ran ge of possible values for Qo 
(which we know from section 
than 10.0). 



II C must have a value less 



In all these data sets, the horizon remains on the throat 
until some critical value of Qo is reached and then it de- 
parts and moves outward in coordinate space. Until the 
horizon leaves the throat, Cp/CA continues to grow, but 
as soon as it departs this ratio begins to decrease. This 
downward trend in Cp/CA continues as Qo is increased 
until a value of Qq is reached which is inconsistent with 
asymptotic flatness. The most distorted spacetime oc- 
curs when the largest value of a is used, although the 
growth is very slow after the value of exp(i7) has increased 
to about half the fundamental £ = 2 wavelength of the 
black hole. 

While the information we obtained does not provide us 
with a definite answer about the maximum possible pro- 
lateness of the apparent horizon in these spacetimes, it 
suggests that there is some limiting value of Cp/CA- The 
largest value we have measured is Cp/CA — 88. However, 
it is possible that despite the observed trends the quan- 
tity Cp/Ca may be made to grow without limit for some 
untried Brill wave shape or combination of shapes. 



H. Maximum Oblateness 

How oblate can a black hole be? To get a feel for the 
meaning of the ratio Cp/CA we note that for an infinitely 
thin disk in Euclidean space we have Cp/CA — 0.900 and 
for an extremal Kerr black hole (a/m = 1) we have a 
smaller value, Cp/CA = 0.860. Thus, although these 
ratios are quite close to unity, they actually represent 
extremely distorted surfaces. For a very wide range of 
Brill wave parameters we have seen a robust lower bound 
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for Cp/Ca ~ 0.77, with Qq = —0.65, as discussed above. 
This seems to follow the trend discussed above for the 
prolate black hole spacetimes. 

I. Penrose Inequality 

In 1973 Penrose proposed a criterion for a spacetime 
which, if violated, would indicate Cosmic Censorship 
would be violated |^^. The criteria is simply that the 
irreducible apparent horizon mass must be less than or 
equal to the ADM mass. 

Jang and Wald proved |^ the inequality for space- 
times with a single AH of topology which also sat- 
isfies a certain mathematical condition. (It essentially 
states that it is possible to foliate the region of the space- 
time exterior to the apparent horizon with a sequence of 
nested 2-surfaces, which are asymptotically spheres at 
large radii, using a "lapse" whose value is equal to the 
inverse of the extrinsic curvature of those 2-surfaces) We 
have not shown that this condition is actually satisfied 
for the wide class of black hole spacetimes considered 
here, but we have seen that this inequality is preserved 
in all spacetimes considered in this paper. Bernstein also 
sought for and failed to find a violation of the Penrose 
inequality 0Jl3l using the Brill wave distortions of the 
non-rotating black hole, but did not consider the extreme 
distortions used in this paper. A violation of the Penrose 
inequality would result in a negative value for the MRL, 
and we can see that this does not occur if we review 
Figs. I and |l|. 

V. SUMMARY AND FUTURE WORK 

In this paper we have described the construction of 
new initial data sets for distorted rotating black holes in 
detail. These data sets correspond to a large family of 
black hole spacetimes, generalizing the Kerr [Q, Bowcn 
and York [|l9|, distorted Schwarzschild black holes 
and should provide a rich testing ground for evolutions 
of highly dynamic, rotating and nonrotating black holes. 

We have analyzed many of the properties of extreme 
parameter choices these data sets. We have, for example, 
located the analytic upper and lower limits to the possi- 
ble distortion parameters for these spacetimes, located 
the transition points beyond which multiple trapped 
surfaces exist, and identified spacetimes with the most 
oblate/prolate apparent horizons in our family of data 
sets. Furthermore, we investigated the effect of rota- 
tion and cosmic-screw type distortions on these data sets. 
This investigation should help provide a roadmap to in- 
teresting, distorted black hole spacetimes for evolution in 
both 2D and 3D. 

In the near future we intend to apply the horizon find- 
ing code described in Ref. |4l| to some of the above data 
sets, to determine to what extent the behavior of the 



event horizons is in agreement with the behavior of the 
more extremely distorted apparent horizons. Do space- 
times with strongly distorted apparent horizons have 
strongly distorted event horizons? This research provides 
us with important data for this investigation. In particu- 
lar, (7 = 2 configurations with amplitudes Qo > 3.5 seem 
likely to yield the most distorted spacetimes. We also 
plan to evolve the cosmic screw spacetimes and study 
them in more detail. 
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